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An analytic formula is given for the traceless transverse part of the anisotropic stress 
tensor due to free streaming neutrinos, and used to derive an integro-differential equa- 
tion for the propagation of cosmological gravitational waves. The solution shows that 
anisotropic stress reduces the squared amplitude by 35.6 % for wavelengths that enter the 
horizon during the radiation-dominated phase, independent of any cosmological param- 
eters. This decreases the tensor temperature and polarization correlation functions for 
these wavelengths by the same amount. The effect is less for wavelengths that enter the 
horizon at later times. At the longest wavelengths the decrease in the tensor correlation 
functions due to neutrino free streaming ranges from 10.7% for S]jm/i 2 = 0.1 to 9.0% for 
^Adh 2 = 0.15. An Appendix gives a general proof that tensor as well as scalar modes sat- 
isfy a conservation law for perturbations outside the horizon, even when the anisotropic 
stress tensor is not negligible. 

I. Introduction 

It is widely expected that the observation of cosmological tensor fluctua- 
tions through measurements of the polarization of the microwave background 
may provide a uniquely valuable check on the validity of simple inflationary 
cosmologies. For instance, for a large class of inflationary theories with single 
scalar fields satisfying the "slow roll" approximation, the wave-number de- 
pendence Vs oc /c™ s_1 and Vt k nT of the scalar and tensor power spectral 
functions and the ratio of these spectral functions after horizon exit during 
inflation are related by[l] 

VrjVs = -n T /2 . (1) 

But in order to use observations to check such relations, we need to know 
what happens to the fluctuations between the time of inflation and the 



present. There is a very large literature on the scalar modes, but ever since 
the first calculations [2] of the production of tensor modes in inflation, with 
only one exception[3] known to me, the interaction of these modes with mat- 
ter and radiation has simply been assumed to be negligible in studies of the 
cosmic microwave background [4]. It is not included in the widely used com- 
puter program of Seljak and Zaldarriaga[5]. As we shall see, the effect is not 
negligible even at the relatively low values of I where the B-type polarization 
multipole coefficients Cbi are likely to be first measured, and becomes quite 
significant for larger values of I. 



II. Damping Effects in the Wave Equation 




The interaction of tensor modes with matter and radiation vanishes in the 
case of perfect fluids, but not in the presence of traceless transverse terms in 
the anisotropic stress tensor. In general, the tensor fluctuation satisfies 

= 167^ , (2) 

where dots indicate ordinary time derivatives. Here the components of the 
perturbed metric are 

0oo = -1 , 9io = , 9ij(x, t) = a 2 {t) + /^-(x, *)] (3) 

where ^ (x, t) is treated as a small perturbation; and ^ (x, t) is the anisotropic 
part of the stress tensor, defined by writing the spatial part of the perturbed 
energy-momentum tensor as = p + a 2 7iij , or equivalently 

M, • (4) 

where p is the unperturbed pressure. In these formulas we are considering 
only tensor perturbations, so that 

ha = , dihij = , iru = , 8^ = . (5) 

For a perfect fluid -k^ = 0, but this is not true in general. For instance, in 
any imperfect fluid with shear viscosity rj, we have [6] = —rjhij. Neverthe- 
less, as we shall show in the Appendix, even where hydrodynamic approxi- 
mations are inapplicable, hij becomes time- independent as the wavelength 
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of a mode leaves the horizon, and remains time-independent until horizon 
re-entry. All modes of cosmological interest are still far outside the horizon 
at the temperature ~ 10 10o K where neutrinos are going out of equilibrium 
with electrons and photons, so can be effected by anisotropic inertia only 
later, when neutrinos are freely streaming. 

We can calculate the contribution of freely streaming neutrinos to 
exactly [7]. We define a density n(x, p,t) as 



n(x, p,*)=£ n - 4(0) n ^ - prM) 



(6) 



\i=i 



\i=l 



with r labeling individual neutrino and antineutrino trajectories. The rela- 
tivistic equations of motion in phase space for any metric with g 00 = — 1 and 
g i0 = are 

p* PiPr ( &9jk\ 



P°r 



Pr 



2p° V dx i 



It follows then that n satisfies a Boltzmann equation 

dn p l dn pip k dg^ 



dn 
~dt 



0, 



(7) 



(8) 



dx l p° dpi 2p° dx l 

it being understood that p % and p° are expressed in terms of the independent 
variable pi by p % = g^Pj and p° = (g^PiPj) 1 ^ 2 - At a time t\ soon after 
neutrinos started free streaming, n had the ideal gas form (assuming zero 
chemical potentials) 



n(x,p,*i) = 



N 



exp UgVfctJpiPj/kBTi) + 1 



= m(x,p), (9) 



where is the number of types of neutrinos, counting antineutrinos sepa- 
rately, and he is Boltzmann's constant. We therefore write 



n(x, p, t) = ni(x, p) + <5n(x, p, t) 



(10) 



so that 5n vanishes for t — t\. 

In the absence of metric perturbations, Eq. (8) and the initial condition 
(9) have the solution n(p) = n(p), where n{p) is the zeroth-order part of n±: 



n(p) 



N 



exp {p/ksTiax) + 1 



(11) 
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and p = y/piPi. To first order in metric perturbations, Eq. (8) gives 
afn(x,p,f) _^ flfa(x,p,t) P—n'( p )fi.fi.fi.JL ( h ..( x t ) - Mx 

(12) 

where hats denote unit vectors. (In putting the Boltzmann equation in this 
form, we use that fact that n\ depends on x and pi only through the combina- 
tion <7 u '(x, ti)piPj, so that to first order <9fc7ii(x, p) = —p n' '(p)pipjdkhij(x., t±).) 

We now suppose that the x-dependence of ^ (x, t) is contained in a factor 
exp(ik • x), where k is a co- moving wave number. 1 Eq. (12) and the initial 
condition that Sn — at t — ti then have the solution 

Sn(p,u) = - l -pn\p)p- kpipj £du' ^ ~^ (h^u') - ^(0)) (13) 

where we now drop the position argument, and write Sn and hij as functions 
of a variable u instead of t, with u defined as the wave number times the 
conformal time 

f* dt' 



u = k —— . (14) 
Jti a(t') y ' 

The space part of the neutrino energy-momentum tensor is given by 
^ = -7^= E %* 3 (* " x,) = ~J= i ( A d Ph ) ^ (15) 



This yields terms of first order in hij(u) from p % = g' l ^Pj and p° = yg^PiPj, a 
term of first order in hij(0) from the term n\ in n, and a term of first order 
in hij(u) — hij(0) from Sn. Collecting all these terms and using Eq. (5) yields 
a surprisingly simple formula for n^: 

7Ti» = -4p„(u) j K(u - U) h'^iU) dU (16) 

where primes now indicate derivatives with respect to U or u; K is the kernel 

tw x 1 f +1 r /-, 2\2 .- TS sins 3coss 3 sins 

K(s) S -J d X (l - *>)V" = - - + (17) 



1 Conventionally the co- moving coordinate x and wave number k are normalized by 
defining a(t) so that a = 1 at present. Here we will leave this normalization arbitrary. 
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and p v = a~ A J d 3 ppn(p) is the unperturbed neutrino energy density. 

To continue, we use Eq. (16) in Eq. (2) and express time-derivatives 
in terms of -u-derivatives. This gives an integro-differential equation for 

hij(u)[8}: 

Kjiu) + 2 -^K J {u) + M«) = -24/„(u) f^f) [ U K{u - U) KjiU) dU , 

(18) 

where /„ = p v /p. 

We took the initial time t\ to be soon after neutrinos started free stream- 
ing, so interesting perturbations are outside the horizon then, and for some 
time after. As we show in the Appendix, hij rapidly became time indepen- 
dent after horizon exit, and remained so until horizon re-entry. In terms of 
u, we then have the initial condition 

K J (0)=0. (19) 
The solution of Eq. (18) can therefore be put in the general form 

hij(u) = MO)x(u) (20) 
where x( u ) satisfies the same integro-differential equation as hij(u) 

+ ^-x'H + X (u) = -24/„(u) (^) [ U K{u - U) X \U) dU; , 
a[u) \ a[u) J Jo 

(21) 

and the initial conditions 

X (0) = 1 , x'(0) = . (22) 



III. Short Wavelengths 

We will first consider wavelengths short enough to have re-entered the 
horizon during the radiation-dominated era (though long after neutrino de- 
coupling), and then turn to the general case in the next section. We can 
take the initial time t\ to be early enough so that it can be approximated as 
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ti ~ 0, with the zero of time defined so that during the radiation-dominated 
era we have a oc \fi. Then in Eq. (21) we can set a' /a = 1/u, while for 3 
neutrino flavors /„ takes the constant value f u (0) = 0.40523. Then Eq. (21) 
becomes 

X » + - X '(u) + x{u) = [ U K{u - U) X '(U) dU , (23) 

u u A Jo 

Because of the decrease of the factor 1/u 2 , the right-hand of Eq. (23) becomes 
negligible for u 3> 1, so deep inside the horizon the solution of Eqs. (22) and 
(23) approaches a homogeneous solution 

X(u) -> A sin(w + 5) ju (24) 

as compared with the solution sin(w)/-u for f u — 0. A numerical solution 
of Eqs. (22) and (23) shows that x{ u ) follows the f v = solution pretty 
accurately until u ~ 1, when the perturbation enters the horizon, and there- 
after rapidly approaches the asymptotic form (24), with A = 0.8026 and 5 
very small. This asymptotic form provides the initial condition for the later 
period when the matter energy density becomes first comparable and then 
greater than that of radiation, so the effect of neutrino damping at these 
later times is still only to reduce the tensor amplitude by the same factor 
A = 0.8026. Hence, for wavelengths that enter the horizon after electron- 
positron annihilation and well before radiation-matter equality, all quadratic 
effects of the tensor modes in the cosmic microwave background, such as 
the tensor contribution to the temperature multipole coefficients Ce and the 
whole of the "B-B" polarization multipole coefficients Cg b, are 35.6% less 
than they would be without the damping due to free- streaming neutrinos. 
(Photons also contribute to 7Tjj, but this effect is much smaller because at 
last scattering photons contribute much less than 40% of the total energy.) 

IV. General Wavelengths 

To deal with perturbations that may enter the horizon after the matter 
energy density has become important, let us switch the independent variable 
from u to y = a(t)/aEQ, where oeq is a(t) at the time £eq of radiation- matter 
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equality. To see how they are related, note that 

^ = —^j- = -^-rHoJsiM (^) 3 + (Oy + n v ) ( fraca a) 4 (25) 
du a EO k a a EO k V V a ) 



The redshift of matter-radiation equality is given by 1 + zeq = ao/ a EQ 
fi M /(fi 7 + fi„), so Eq. (25) can be simplified to read 



du Q 



dy V l + V 

where 



(26) 



Q = k . (27) 

a H ^JVL M (l + z E q) 



Since u — >■ for y — > 0, the solution of Eq. (26) is 



?! 



2Q( v /l+^-l). (28) 

The Hubble constant at matter-radiation equality has the value H E q = 
H ^2Q M {^ + z eq) 3 , so Eq. (27) can be written 

Q = V2k/k EQ , (29) 

where k E Q = cl E qH E q is the wave number of perturbations that just enter 
the horizon at the time of radiation-matter equality. (Hence in particular 
the results of the previous section apply for Q 1.) 

The fraction of the total energy density in neutrinos is well known 

, / a n u (a /a) A = /„(0) 

Myj fi M (ao/a) 3 + (fi 7 + Ome^)(a /fl) 4 1 + y 1 J 

where 

^n^T ' 40523 ' (31) 

A little algebra then lets us put Eq. (21) in the form 

(32) 
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where K(y,y') is the same as the K(s) given by Eq. (17), but with s now 
given by 

s = z -z> = 2Q( V /TT^ - y/l+y') (33) 
The initial conditions (22) now read 



x(0) ^ dxiv) 



dy 



= . (34) 



We now have to face the complication that for general Q the value of y at 
last scattering is not in an asymptotic region where the effect of anisotropic 
inertia is simply to damp x(t) by some constant factor. We therefore now 
have to consider what feature of x(t) is related to observations of the cosmic 
microwave background. It is x that enters into the Boltzmann equation 
for perturbations to the temperature and Stokes parameters [9], so in the 
approximation of a sudden transition from opacity to transparency, we expect 
the B-B and other multipole coefficients to depend on x(v) only through a 
factor \x'{Ul)\ 2 , where y^ — (1 + zeq)/(1 + zl) is the value of y at last 
scattering. Hence we will be primarily interested in calculating the value of 
\x\Ul)\ 2 for various values of Q, and comparing these values with what they 
would be in the absence of anisotropic inertia. 

For T 70 = 2.725°K, we have fl^ + tt u = 4.15 x lO" 5 /*" 2 , so, taking l + z L = 
1090, the parameter y L is 

y L = 22.m M h 2 . 

It will be useful also to have an idea of the value of £ for which the multipole 
coefficients in various correlation functions are dominated by perturbations 
with a given Q. The dominant contribution to a multipole coefficient of 
order £ comes from wave numbers k ~ a-L^/di,, where is a(t) at the time 
of last scattering, and di is the angular diameter distance of the surface of 
last scattering, which for flat geometries is: 

dx 

a ; 



— /' 

- + z T ) hi 



H (l + z L ) Ji/(i+z L ) ^q mX + (i _ Q M ) x i ' 



where zl is the redshift of last scattering. Thus the multipole order that 
receives its main contribution from wave lengths that are just coming into 
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the horizon at matter-radiation equality is 



d L k E Q rzTL — — r I' 1 dx 



£ EQ = = J2Q M (1 + z m ) / , (35) 

O.L V ^ h/(l +ZL ) + (1 _ fi M ) x 4 

where zeq is the redshift of matter-radiation equality. For present radiation 
temperature T 70 = 2.725°K and f2 M ^ 2 = 0.15 this redshift is z E q = 3613. 
If also Vl M = 0.3 and 1 + z L — 1090 then the integral in Eq. (35) has the 
value 3.195, and so Eq. (35) gives £ E q = 149. Hence for these cosmological 
parameters, Eq. (29) gives 

P P 

Q = ^j- - 4 • 

«EQ 105 

When referring below to specific values of I, it will always be understood 
that the conversion from Q to I has been made using these cosmological 
parameters, but it should be kept in mind that the dependence of the function 
x(y) on y and Q is independent of cosmological parameters, and that the 
value of y at last scattering depends only on T 70 , 1 + zl- and flMh 2 , not on 
VLm or f2 vac . 

Let us first consider the case Q <^ 1, which for the above cosmological 
parameters corresponds to £ 100. Here the kernel K(y, y') has the constant 
value 1/15, and Eqs. (32) and (34) have a solution of the form 

x (y) - 1 - Q 2 g(y) for Q - (36) 

where g(y) is independent of Q, and satisfies the inhomogeneous differential 
equation 

,d 2 g(y) {2(l + y) l\ dg(y) 8f„(0) , s , 



and the initial conditions 



^(0) = ^'(0) = . (38) 



According to the above discussion, the streaming of free neutrinos damps 
the various tensor correlation functions of the cosmic microwave background 
by a factor \x' {u l) / Xo{y l)\ 2 , which for Q < 1 becomes \g' {y l) / g'^y l)\ 2 , the 
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subscript denoting quantities calculated ignoring this damping, i.e., for 
/„ = 0, and ul again equal to the ratio of a(t) at last scattering to that 
at matter-radiation equality. Numerical solutions of Eqs. (37) and (38) for 
/i/(0) = 0.40523 and for f u = show that the damping factor \g'(yL)/ 9o(yL,)\ 2 
is very close to a linear function of yi and hence of I^m^ 2 for observationally 
favored values of ^m^ 2 , increasing from 0.893 at f2^/i 2 = 0.10 to 0.910 for 

n M h 2 = o.i5. 

This damping is relatively insensitive to Q for small Q. For instance, 
numerical integration of Eqs. (32) and (34) shows that for f2^-/i 2 = 0.15, the 
damping has only decreased from 9% to 8% for Q = 0.55 (£ ~ 58), and to 7% 
for Q = 0.8 (£ ~ 84). Matters are more complicated for larger values of Q and 
£, because the damping factor |x'(|/l)/Xo(|/l)| 2 is the ratio of two oscillating 
functions with slightly different phases, so that the plot of \x' {y l) / x'oiy l)\ 2 
vs. Q shows narrow spikes: this ratio becomes infinite at values of Q for 
which Xo(Ul) vanishes and then almost immediately drops to zero at the 
slightly larger value of Q for which x'{yL) vanishes. (Even if we average 
over the small range of y values over which last scattering occurs, the plot of 
{W {y l) / 'x'oiy l)\ 2 ) vs. Q still shows finite though high narrow spikes at the 
zeroes of x'o(i/l)-) These spikes are not particularly interesting, because they 
occur at values of Q where x'{Ul) is particularly small, so that the multipole 
coefficients in the various tensor temperature and polarization correlation 
functions will be very difficult to measure for the corresponding values of 
i. The values of \x' {y l) / x'oiy l)\ 2 in the relatively flat regions between the 
spikes steadily decreases from the value ~ 0.9 for Q 1 to a value close to 
the result .644 found in the previous section for Q — 10. 

The effects considered in this paper will doubtless eventually be taken 
into account in the computer programs used to analyze data from PLANCK 
and other future facilities. In the meanwhile, the planning of future observa- 
tions should take into account that the damping of cosmological gravitational 
waves is not negligible. 
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APPENDIX: SUPERHORIZON CONSERVATION LAWS 

This Appendix will prove a result quoted in Section II, that in all cases 
there is a tensor mode whose amplitude remains constant outside the hori- 
zon, even where some particles may have mean free times comparable to 
the Hubble time. The argument is similar to one used previously to show 
the existence under very general conditions of two scalar modes for which a 
quantity 1Z is constant outside the horizon. [10] It is based on the observa- 
tion that for zero wave number the Newtonian gauge field equations and the 
dynamical equations for matter and radiation as well as the condition k = 
are invariant under coordinate transformations that are not symmetries of 
the unperturbed metric. 2 The most general such transformations are 

x° ->ar° + e(*);, x i -> (Sij - -u^ x j , (Al) 

where H = a/a, e(t) is an arbitrary function of time, and lo^ = ujji is 
an arbitrary constant matrix. Under these conditions we have something 
like a Goldstone theorem: since the metric satisfies the field equations both 
before and after the transformation, the change in the metric under these 
transformations must also satisfy the field equations. This change is simply 

8 goo = e(f) , Sg i0 = , Sg i:j = a 2 (t) 

This means that for zero wave number we always have a solution with scalar 
modes 

= He - tUii/3 , $ = -e (A3) 

and a tensor mode 

= uj^ - -5ijUj kk . (A4) 

2 In this respect, the theorem proved here is similar to the Goldstone theorem[ll] of 
quantum field theory. The modes for which 1Z or hij are constant outside the horizon take 
the place here of the Goldstone bosons that become free particles for long wavelength. 



- H{t)e(t)8ij + uj. 



(A2) 
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(The notation for $ and \I/ is standard, and the same as in Ref. [10].) These 
are just gauge modes for zero wave number, but if they can be extended to 
non-zero wave number they become physical modes, since the transforma- 
tions (Al) are not symmetries of the field equations except for zero wave 
number. For the scalar modes there are field equations that disappear in the 
limit of zero wave number, so that conditions $ = ^ — 8nGirs and Su = e 
(where ns is the scalar part of the anisotropic inertia, called a in Ref. [10]) 
and Su is the perturbation to the velocity potential) must be imposed on the 
solutions (A3) for them to have an extension to non-zero wave number. It 
follows then that the zero wave number scalar modes that become physical 
for non-zero wave number satisfy 

e = —He + — StcGtts , Su = e . (A5) 

Then for zero wave number the quantity 1Z = — \I/ + hSu has the time- 
independent value 

K = u kk /3; . (A6) 

For tensor modes there are no field equations that disappear for zero wave 
number, so the solution hij =constant automatically has an extension to a 
physical mode for non-zero wave number. 

As examples, we note that both the anisotropic stress tensor 71"^ = —rjhij 
for an imperfect fluid with shear viscosity 77 and the tensor (16) for freely 
streaming neutrinos vanish for = 0, so in the limit of zero wave numbers 
Eq. (2) has a solution with = 0. The above theorem shows that this result 
applies even when some particle's mean free time is comparable with the 
Hubble time, in which case neither the hydrodynamic nor the free- streaming 
approximations are applicable. 

The solution with = for zero wave number is not the only solu- 
tion, but the other solutions decay rapidly after horizon exit. There is no 
anisotropic inertia in scalar field theories, and in the absence of anisotropic 
inertia, Eq. (2) for zero wave number has two solutions, one with con- 
stant, and the other with oc a -3 , for which rapidly becomes constant. 
The energy-momentum tensor of the universe departs from the perfect fluid 
form later, during neutrino decoupling, and perhaps also during reheating 
or periods of baryon or lepton nonconservation, but during all these epochs 
cosmologically interesting tensor fluctuations are far outside the horizon, and 
hence remain constant. 
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